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Abstract 

This  report  gives  a  review  of  the  various  types  of  unsteady  flow  phenomena  that  occur  within 
a  turbomachine  and  examines  existing  methods  for  computing  several  of  these  flows.  Then, 
we  discuss  a  method  for  applying  two-dimensional,  unsteady,  thin  airfoil  theory  in  a  stripwise 
manner  to  compute  the  unsteady  forces  and  moments  on  the  shaft  of  a  rotor  operating  within 
a  nonuniform  inflow.  Techniques  for  computing  these  unsteady  forces  and  moments  on  a  rotor 
blade  row  with  rake  and  skew  are  included. 
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Introduction 

For  many  years,  engineers  have  designed  and  analyzed  turbomachines  using  one  major  as¬ 
sumption.  They  have  assumed  the  flow  to  be  steady.  However,  as  we  strive  to  design  even  better 
turbomachines,  we  must  begin  to  focus  on  the  complexities  of  unsteady  flow.  For  various  types 
of  turbomachinery,  unsteady  flow  can  affect  energy  transfer,  blade  vibration,  noise  generation, 
heat  transfer,  and  cavitation  characteristics. 

Before  we  discuss  the  unsteady  flow  considerations  in  turbomachinery  design,  we  will  briefly 
review  the  various  types  of  unsteady  flows  in  turbomachines.  Given  some  knowledge  of  the 
physical  phenomena  involved,  we  shall  examine  several  available  methods  for  calculating  unsteady 
forces  and  moments  and  how  these  methods  can  be  applied  to  turbomachinery  blading.  We 
shall  discuss  the  advantages  and  disadvantages  of  all  these  methods  as  well  as  point  out  the 
approximations  needed  in  the  development  of  these  methods.  Finally,  we  show  how  we  can 
adapt  these  methods  to  turbomachinery. 
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Unsteady  Flows  in  Turbomachines 


Before  we  include  unsteady  effects  in  our  design  process,  we  must  begin  to  understand  the 
physics  of  these  various  unsteady  effects  and  the  relationship  between  these  various  effects.  Dean 
[1959]  showed  that  the  flow  through  a  turbomachine  had  to  be  unsteady.  His  “unsteady  paradox” 
states  that  local  unsteady  changes  in  static  pressure  must  exist  within  a  turbomachine  in  order  for 
that  machine  to  exchange  energy  with  the  fluid.  If  this  turbomachine  includes  only  a  single  blade 
row  with  a  uniform  inflow,  we  can  consider  the  problem  to  be  steady  by  writing  our  governing 
equations  in  a  coordinate  system  rotating  with  the  blades.  When  a  row  of  rotating  blades 
operates  with  a  circumferentially  varying  inflow,  or  if  it  operates  in  the  vicinity  of  stationary 
blades,  struts,  or  other  appendages-the  flow  is  unsteady,  despite  the  frame  of  reference  we  use. 

When  we  discuss  unsteady  flow  in  turbomachinery,  we  usually  refer  to  periodic  unsteadiness. 
We  must  also  remember  that  a  random  or  aperiodic  unsteadiness  exists  because  of  turbulence. 
Turbulence  also  has  strong  effects  on  energy  and  heat  transfer  and  it  also  includes  important 
random  pressure  fluctuations.  Our  concern  with  turbulence  will  be  in  how  turbulence  interacts 
with  periodic  unsteadiness. 

Let  us  now  discuss  why  an  inherent  periodic  unsteadiness  exists  within  a  turbomachine. 
Blades  moving  relative  to  one  another  can  operate  within  a  periodically  varying  potential  flow 
field.  For  incompressible  flows,  we  can  normally  neglect  this  periodic  potential  interaction  when 
a  blade  row  is  displaced  at  least  one  chord  length  upstream  or  downstream  of  neighboring  rows 
of  blades  or  appendages.  For  compressible  flows,  the  influence  of  the  periodic  potential  field 
extends  much  farther  axially  as  the  Mach  number  increases. 

Because  of  the  fairly  large  axial  blade  gaps  in  most  turbomachines,  the  periodic  potential 
interactions  are  usually  dominated  by  other  types  of  periodic  unsteadiness.  Most  importantly, 
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a  blade  row  interacts  with  the  wakes  shed  from  an  upstream  row  of  blades  or  appendages  that 
have  relative  motion.  These  wakes  are  convected  downstream,  chopped  by  the  downstream 
blade  row,  and  transported  through  the  downstream  blade  passages.  As  seen  in  Figure  1,  the 
periodically  passing  wakes  cause  a  periodically  varying  incidence  angle  on  the  downstream  blades. 
This  unsteady  incidence  results  in  a  periodically  varying  blade  static-pressure  distribution  and 
thus  a  periodically  varying  lift.  In  addition,  the  behavior  of  the  chopped  wake  segment  as  it 
is  transported  through  the  blade  passage  is  quite  interesting.  For  instance,  we  can  see  from 
Figure  2  that  the  highly  turbulent  fluid  in  the  wake  segment  moves  from  the  suction  surface  of  a 
compressor  blade  and  impinges  on  the  pressure  surface  as  the  wake  segment  moves  downstream. 
This  action  results  in  local  pressure  fluctuations  due  to  the  passing  wake  segment.  Meyer  [1958] 
made  the  first  theory  of  this  wake  interaction  effect  while  Kerrebrock  and  Mikolajczak  [1970] 
proved  it  experimentally.  For  a  turbine  blade  passage,  the  wake  fluid  moves  from  the  pressure 
surface  to  the  suction  surface.  Other  researchers  have  further  extended  our  knowledge  of  wake 
interaction  effects. 

Periodic  wake  interactions  have  many  influences  that  can  affect  a  turbomachine  design.  The 
periodically  varying  blade  static-pressure  distribution  can  cause  a  forced  vibration.  Unsteady 
forces  and  moments  affect  the  blades  and  whatever  is  attached  to  the  blades  such  as  the  rotor 
shaft  or  the  stationary  hub  and  casing.  These  unsteady  forces  can  cause  vibration  at  the  frequency 
of  unsteadiness-namely  the  blade  rate  frequency  and  its  harmonics.  Periodic  wake  interactions 
also  cause  the  boundary  layers  near  the  blade  surface  to  be  unsteady.  Varying  the  blade  static- 
pressure  distribution  changes  the  local  streamwise  pressure  gradients  which  in  turn  changes 
the  boundary  layer  characteristics.  In  addition,  the  effects  of  the  transported  wake  segments 
can  have  strong  local  effects  on  the  boundary  layer.  Evans  [1978],  Walker  [1982],  and  Hodson 
[1983]  all  made  measurements  within  turbomachines  that  show  some  boundary  layers  alternating 
between  laminar  and  turbulent  because  of  the  unsteady  flow.  Periodic  unsteadiness  can  also  affect 
separation.  These  unsteady  boundary  layers,  along  with  the  complexities  of  the  wake  mi  eractions, 
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strongly  affect  the  energy  transfer  within  the  machine.  In  addition,  the  changing  character  of 
the  boundary  layers  affects  the  turbulence  and  thus  the  random  pressure  fluctuations.  Finally,  if 
the  pressure  fluctuations  become  large  enough,  cavitation  can  occur  in  machines  operating  with 
liquid  as  the  working  fluid-as  shown  by  the  unsteady  cavitation  experiment  of  Yeh  and  Eisenhuth 
[1959}. 

A  few  other  characteristics  of  periodic  wake  interactions  should  be  noted.  First,  time  or 
phase  lags  always  appear.  As  the  incidence  angle  changes,  the  static  pressure  near  the  leading 
edge  of  the  blade  changes  much  sooner  than  the  static  pressure  near  the  trailing  edge.  In  fact, 
experiments  have  shown  that  the  amplitude  of  the  periodically  changing  pressure  is  much  larger 
near  the  leading  edge  than  near  the  trailing  edge.  The  effect  of  the  varying  incidence  appears 
to  decay  as  one  approaches  the  trailing  edge.  In  addition,  the  characteristics  of  the  unsteady 
boundary  layers  can  also  lag  behind  the  change  in  the  local  velocity  and  static-pressure  gradient. 
We  should  also  take  note  of  a  second  important  factor.  The  flow  within  the  wakes  and  wake 
segments  can  also  be  unsteady.  We  have  all  seen  measurements  that  make  us  think  of  a  wake 
as  merely  a  region  of  velocity  defect.  However,  these  measurements  are  probably  time-average 
measurements.  In  fact,  real  turbomachine  wakes  can  be  made  of  periodically  unsteady  vortex 
streets  such  as  those  we  normally  identify  with  the  wakes  behind  blunt  bodies.  This  wake 
unsteadiness  includes  time  varying  pressures  which  can  influence  the  unsteady  wake  interactions 

In  addition  to  the  periodically  unsteady  potential  and  wake  interactions,  we  can  encounter 
other  kinds  of  periodically  unsteady  flows.  Blades  may  also  periodically  interact  with  shock 
waves,  end-wall  boundary  layers,  tip  vorticies,  and  other  secondary  flow  vorticies.  The  three- 
dimensional  nature  of  these  types  of  flow  can  create  significant  circumferential  variations  in 
velocity  and  static  pressure  which  interact  unsteadily  with  downstream  blades. 

We  have  only  discussed  unsteady  flows  that  are  inherently  found  in  turbomachine  flows. 
Since  these  flows  will  always  exist,  we  must  learn  their  effects  and  use  this  knowledge  to  design 
better  turbomachines.  Other  types  of  unsteady  flows  can  occur,  however,  which  we  must  try  to 
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avoid.  Thes'  unsteady  flows  are  instabilities  which  may  cause  large  resonant  stresses  in  blades 
and  drastic  reductions  in  machine  performance. 

As  blades  vibrate,  the  mechanical  surfaces  exchange  energy  with  the  fluid  passing  through 
the  blade  passage.  If  the  fluid  extracts  energy  from  the  blade,  the  surface  oscillations  will  be 
aerodynamically  dampened.  However,  if  the  fluid  increases  the  energy  of  the  oscillating  blades, 
the  blades  will  become  vibrationally  unstable.  We  call  this  instability  flutter.  Flutter  exists 
because  the  unsteady  work  on  the  blades  exceeds  the  work  dissipated  by  friction  and  material 
damping.  Normally,  flutter  occurs  when  the  blade  boundary  layer  separates.  But  flutter  can 
also  occur  when  the  blade  boundary  layers  are  attached-especially  at  near-choke  conditions. 
Compressibility  strongly  increases  the  chances  of  developing  flutter. 

Large  inlet  flow  distortions  can  cause  a  blade  to  stall  in  a  turbomachine  in  a  similar  manner 
in  which  an  isolated  blade  may  stall.  In  both  cases,  the  potential  flow  around  the  blade  will 
change.  In  a  turbomachine,  however,  the  new  potential  flow  will  decrease  the  incidence  of  one 
neighboring  blade  and  will  increase  the  incidence  of  the  other  adjacent  blade,  causing  this  blade 
to  stall.  A  stall  zone  of  one  or  more  blade  passages  will  then  propagate  circumferentially  around 
the  annulus  of  blades.  We  call  this  flow  instability  rotating  stall  and  it  leads  to  a  large  pressure 
drop  across  the  blade  passage.  Hysteresis  effects  cause  the  recovery  from  rotating  stall  to  be 
very  difficult.  For  machines  where  the  working  fluid  is  liquid,  inlet  flow  distortions  and  the 
corresponding  changes  in  incidence  angle  can  lead  to  cavitation.  Some  machines  have  shown  an 
instability  similar  to  rotating  st  ill  called  rotating  cavitation. 

Rotating  stall  is  one  factor  that  can  lead  to  a  global  instability  within  a  compressor  called 
surge.  Surge  occurs  when  the  throttle  pressure-loss  characteristic  curve  has  a  smaller  slope  than 
the  constant-speed  compressor  characteristic  curve  at  the  point  of  operation.  A  reduction  in 
flow  would  then  cause  a  greater  reduction  in  compressor  pressure  rise  than  the  corresponding 
pressure  drop  through  the  throttling  system.  As  a  result,  the  annulus  mass  flow  oscillates  with 
large  changes  in  amplitude  at  low  frequencies. 
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Reduced  Frequency 


Periodic  unsteady  flow  can  have  strong  effects  on  a  turbomachine,  but  the  strength  of  these 
effects  varies  with  the  level  of  unsteadiness.  Therefore,  we  need  some  parameter  that  represents 
the  level  of  unsteadiness. 

The  level  of  unsteadiness  depends  on  two  time  scales.  First,  the  time  scale  associated  with 
the  angular  frequency  of  the  unsteady  motion  relative  to  the  rotating  blades  is  •£-.  For  this 
case,  our  interest  is  in  the  nth  harmonic  of  the  unsteady  motion.  Second,  we  need  the  time  scale 
associated  with  the  time  it  takes  a  fluid  particle  to  travel  through  the  device.  If  this  device  is  a 
blade  row,  a  fluid  particle  will  travel  through  the  blade  passage  in  time  ~.  The  ratio  of  these 
time  scales  then  describes  the  change  in  local  flow  parameters  while  a  fluid  particle  resides  in  a 
blade  passage.  If  we  change  the  relevant  length  from  c  to  |,  we  refer  to  this  ratio  as  the  reduced 
frequency, 

_  vnc 
W"  ~  2  W  ' 

Sometimes  c  is  used  instead  of  We  could  also  think  of  the  reduced  frequency  as  a  ratio  of 
length  scales-the  semi-chord  divided  by  the  wavelength  of  the  n,h  harmonic  of  the  unsteady 
motion. 

For  many  unsteady  flow  problems,  such  as  the  periodic  shedding  of  vorticies  into  a  wake, 
one  refers  to  the  reduced  frequency  as  a  Strouhal  number  (see  Strouhal  [1878]).  For  the  prob¬ 
lem  of  a  shedding  Karman  vortex  street,  we  might  substitute  the  shedding  frequency  for  un  and 
the  trailing  edge  diameter  for  |  in  defining  a  Strouhal  number.  For  periodic  interaction  prob¬ 
lems  in  turbomachinery,  we  use  the  reduced  frequency  where  the  disturbance  is  both  moving 
circumferentially  relative  to  a  blade  row  and  is  being  convected  downstream. 
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For  a  reduced  frequency  much  less  than  one,  we  refer  to  the  unsteady  flow  as  quasi-steady. 
That  is,  the  flow  is  only  a  function  of  the  instantaneous  conditions  and  not  a  function  of  pre¬ 
vious  conditions.  For  a  reduced  frequency  greater  than  one,  the  limiting  state  of  quasi-steady 
flow  would  provide  estimates  of  the  unsteady  flow  that  would  be  much  too  conservative.  Most 
turbomachinery  flows  are  characterized  by  reduced  frequencies  of  order  unity  or  greater  and  we 
cannot  adequately  describe  these  flows  in  a  quasi-steady  manner.  In  compressible  flow,  the  level 
of  unsteadiness  depends  on  the  product  of  the  reduced  frequency  and  the  Mach  number. 

In  turbomachinery  problems,  we  may  require  an  additional  nondimensional  frequency  param¬ 
eter.  Once  again,  we  need  the  time  scale  associated  with  the  angular  frequency  of  the  unsteady 
motion  relative  to  the  rotating  blades,  Next,  we  need  the  time  scale  associated  with  the 
motion  of  the  blades  relative  to  the  distortion,  Taking  the  ratio  gives 


We  call  this  nondimensional  parameter  the  inter-blade  frequency  or  the  inter-blade  phase  angle. 
This  parameter  accounts  for  the  variation  in  unsteady  interactions  from  blade  to  blade.  Just  as 
we  think  of  w„  as  a  ratio  of  semi-chord  to  disturbance  wavelength,  we  can  think  of  r„  as  a  ratio 
of  blade  spacing  to  disturbance  wavelength. 
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Unsteady  Singularity  Solutions 


Several  different  types  of  methods  exist  for  computing  the  unsteady  forces  and  moments  on 
a  turbomachine  blade.  As  we  have  already  discussed,  we  can  use  quasi-steady  methods  for  flows 
with  a  reduced  frequency  much  less  than  one.  We  assume  that  the  blades  operate  continually 
at  the  instantaneous  flow  conditions.  Then  we  can  calculate  the  unsteady  forces  by  actually 
calculating  the  steady  forces  at  each  instant  of  time.  An  example  is  a  quasi-steady  actuator 
disc  analysis.  One  of  the  oldest  methods  for  computing  unsteady  forces  and  moments  involves 
the  modeling  of  the  blades  with  singularities  and  letting  the  singularities  induce  a  potential 
flow  around  the  blade.  Another  method  involves  the  treatment  of  the  blades  as  boundaries  and 
modeling  the  flow  with  the  governing  partial  differential  equations  of  fluid  mechanics.  These 
partial  differential  equations  are  further  approximated  and  solved  numerically. 

Interest  in  the  use  of  singularity  solutions  for  unsteady  flow  problems  began  in  the  1930’s. 
Theodorsen  [1935]  derived  the  exact  solution  for  a  flat  plate  in  sinusoidally  oscillatory  motion. 
Von  Karman  and  Sears  [1938]  used  the  ideas  of  Theodorsen  [1935]  for  their  unsteady  thin  airfoil 
theory  and  finally  Sears  [1941]  expanded  these  ideas  to  derive  his  unsteady  lift  function.  The 

I 

unsteady  thin  airfoil  theory  of  Sears  [1941]  has  many  similarities  with  steady  thin  airfoil  theory 
for  incompressible  flow.  Sears  [1941]  represented  an  airfoil  as  a  flat  plate  of  zero  thickness  and 
modeled  the  airfoil  with  a  continuous  vortex  distribution  along  its  chord.  Each  vortex  induces  a 
potential  flow  at  every  point  in  the  flow  field  which  we  can  superimpose  upon  the  potential  flow 
from  the  steady  inflow.  Unfortunately,  the  system  of  equations  is  not  closed  since  the  downstream 
conditions  remain  unknown.  We  can  complete  the  solution  and  find  a  unique  lift  by  applying 
the  Kutta  condition;  no  velocity  discontinuity  can  exist  at  the  trailing  edge.  Now  we  adjust  the 
strength  of  the  vortex  distribution  such  that  no  fluid  flows  through  the  airfoil  giving  the  airfoil 
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a  unique  bound  vorticity  and  lift.  For  his  unsteady  problem,  Sears  [1941]  treated  the  inflow 
as  a  gust  being  convected  with  a  velocity,  W.  Figure  3  shows  this  gust  which  we  assume  as  a 
harmonic  function  in  time.  Sears  [1941]  assumed  that  the  amplitude  of  the  disturbance  velocity 
normal  to  the  chord  ,v,  is  small  compared  to  W. 

Sears  [1941]  found  it  necessary  to  expand  on  the  Kutta  condition  when  applying  it  to  un¬ 
steady  flow.  Since  the  inflow  varies  in  time,  the  strength  of  the  vortex  distribution  necessary  to 
satisfy  the  blade  boundary  condition  will  also  vary  with  time;  this  then  causes  the  lift  to  vary 
with  time.  Thus,  a  trailing  vorticity  must  be  shed  into  the  wake.  The  trailing  vorticity  affects 
the  flow  throughout  the  flow  field  which  means  that  the  flow  field  depends  on  the  history  of  the 
periodically  unsteady  flow.  Sears  [1941]  modeled  the  trailing  vorticity  with  a  vortex  distribution 
being  convected  away  from  the  airfoil  in  the  direction  of  the  chord  with  a  velocity,  W. 

Sears  [1941]  developed  an  equation  for  the  unsteady  lift  per  unit  span,  Ln,  for  a  transverse 
disturbance  velocity,  vn,  and  a  reduced  frequency,  u„, 

ln(t)  =  l^±pcWvnS(  W„) 

2  da 

where 

S(w„)  =  C(w„)[70(w„)  -  iJi(uj„)}  +  iJi(w„) 

is  called  the  Sears  Function.  Theodorsen  [1935]  developed 

\  _  _ ~  t  Yl(u>n) _ 

1  n>  ~  [7  l(w„)  -  .FlM  +  i[J o(w„)  -  iYo(«n)] 

which  we  call  the  Theodorsen  Function.  We  calculate  unsteady  lift  functions  of  this  kind  from 
Bessel  functions  using  the  appropriate  reduced  frequency.  We  should  note  that  several  of  the 
singularity  techniques  under  discussion  were  developed  for  isolated,  thin  airfoils  with  infinite 
span  where  the  theoretical  value  of  equals  2n.  Since  we  are  concerned  with  nonisolated 
turbomachinery  blades  with  finite  span  and  thickness,  we  will  let  remain  a  variable.  Osborne 
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[1973]  made  a  further  modification.  He  developed  an  unsteady  thin  airfoil  theory  for  transverse 
gusts  that  includes  compressibility  effects. 


Since  the  inflow  gust  varies  harmonically  in  time,  all  the  time  varying  quantities  will  be 
harmonic  functions  which  we  associate  with  an  amplitude  and  a  phase  angle.  We  can  then  plot 
the  Sears  Function  on  the  complex  plane  as  done  in  Figure  4.  For  a  very  small  wn,  S(u>n)  lies 
near  the  real  axis  giving  a  small  phase  angle  and  showing  why  a  quasi-steady  approach  works. 
For  increasing  un,  the  entire  airfoil  cannot  respond  immediately  to  a  changing  incidence  which 
results  in  a  significant  phase  angle. 

Horlock  [1968a]  extended  the  unsteady  thin  airfoil  theory  to  include  a  steady  angle  of  attack 
and  thus  a  chord  wise  gust,  u„, 

£n(0  =  ^~pcW[v„S(ujn)  -I-  au„T(w„)] 

where 

Tn(w)  =  Jo(w„)  4-  *Jl(wn)  +  S(u>r») 


is  termed  the  Horlock  Function.  Holmes  (see  Horlock  [1968b])  and  Naumann  and  Yeh  [1973] 
further  extended  the  unsteady  thin  airfoil  theory  to  include  camber, 


Ut)  =  \^w 


vnS(w„)  +  aunFa(wn)  +  —  Ff(wn) 


where 

Fa  =  J0(w„)  + 


and 

Ff  =  3  J-l(un)  +  Jo(u>n)  +  C(w„)[J0(wn)  -  -  *2Ji(w„)]  . 


They  obtained  this  equation  by  placing  the  vortex  distribution  on  the  camber  line,  but  with 
a  small  value  of  camber  such  that  the  blade  boundary  condition  could  be  satisfied  along  the 
chord  line.  We  should  note  that  these  new  unsteady  lift  functions  are  merely  added  to  the  Sears 
Function.  To  obtain  this  result,  we  must  assume  small  values  of  steady  lift,  steady  angle  of 


11 


attack,  and  maximum  camber.  We  also  need  to  assume  small  values  of  un  and  vn  as  compared 
to  W.  These  assumptions  allow  us  to  linearize  the  various  potential  flow  equations  and  obtain 
Ln  as  the  summation  of  various  functions. 

Before  we  use  these  versions  of  the  unsteady  thin  airfoil  theory  for  turbomachinery  problems, 
we  must  determine  the  correct  reduced  frequency  associated  with  the  various  Bessel  functions. 
For  the  case  of  a  blade  rotating  through  a  nonuniform  inflow,  Figure  5  shows  the  nth  harmonic 
of  an  inflow  distortion  convected  with  an  absolute  flow  angle  0.  For  0  =  0,  the  distortion  lies 
in  the  circumferential  direction  and  has  a  wavelength  of  (2st ) .  For  a  nonzero  value  of  0,  the 
wavelength  becomes  (~")  sec(/?).  However,  we  must  model  this  distortion  as  a  transverse  gust 
and  a  chordwise  gust  in  order  to  apply  unsteady  thin  airfoil  theory. 

As  shown  in  Figure  5,  we  can  decompose  this  distortion  into  two  components-both  perpen¬ 
dicular  and  parallel  to  the  blade  chord.  In  order  to  determine  the  correct  reduced  frequency 
associated  with  these  two  gusts,  we  first  must  compute  the  angular  frequency, 

phase  velocity 
""  =  ~w»v.kngtb  • 

The  phase  velocity  for  the  transverse  (or  perpendicular)  gust  is  the  component  of  the  relative 
velocity  vector,  W,  in  the  chordwise  direction.  As  mentioned  earlier,  the  chorwise  (or  parallel) 
gust  will  only  produce  an  unsteady  lift  if  the  blade  experiences  a  steady  angle  of  attack  or  if 
the  blade  has  camber.  The  phase  velocity  will  be  the  component  of  W  perpendicular  to  the 
chord.  When  a  is  zero,  this  will  give  a  phase  velocity  of  zero  and,  therefore,  values  of  zero  for  the 
angular  frequency,  reduced  frequency,  and  unsteady  lift-which  we  know  cannot  be  the  case  for  a 
cambered  blade.  To  circumvent  this  problem,  we  can  determine  an  effective  angle  of  attack,  aeff, 
by  adding  the  angle  of  zero  lift  to  the  angle  of  attack.  The  zero  lift  line  can  be  approximated  as 
passing  through  the  trailing  edge  and  the  camber  line  at  midchord.  Further,  assuming  that  the 
value  of  camber  at  midchord  is  ym  yields 

o>fr  =  q  +  tan  I  1  . 
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We  can  compute  a  wavelength  for  the  transverse  gust  of  (^-)  csc(7)  and  the  chordwise  gust 
of  - )  sec(7).  FVom  this  information,  we  can  now  find  the  angular  frequency  of  the  unsteady 
transverse  gust  relative  to  the  rotating  blade, 


Vn  =  27T 

and  a  reduced  freqency, 


W  cos(aeff) 


(*?)  csc(7) 


=  sin(7)  cos(aeff)  , 


«n  =  (7^)  sin(7)  cos(aeff)  . 


The  singular  frequency  of  the  unsteady  chordwise  gust  relative  to  the  rotating  blade  is 


i/„  -  2  ir 

and  the  reduced  freqency  is 


W  sin(oeff)  ]  / nW\  ,  x  .  „ 


(  Cos(r)  sin(aefr) 


In  evaluating  the  various  lift  functions,  we  need  to  use  the  appropriate  reduced  frequency  corre¬ 
sponding  to  either  the  transverse  or  chordwise  gust. 


Aside  from  the  assumptions  already  stated,  Horlock  [1968a],  Holmes  (see  Horlock  [1968b]), 
and  Naumann  and  Yeh  [1973]  also  assumed  that  the  incoming  gust  would  not  change  as  it  was 
convected  past  the  airfoil.  Goldstein  and  Atassi  [1976]  allowed  the  gust  to  be  distorted  by  the 
steady-state  potential  flow  field  about  the  airfoil.  This  coupling  effect  led  to  a  second-order 
expansion  for  calculating  the  unsteady  lift  about  the  airfoil.  Atassi  [1984]  expanded  these  ideas 
to  show  that  the  resulting  expressions  for  the  unsteady  lift  still  possess  a  linear  character. 


We  have  discussed  several  methods  to  model  the  unsteady  wake  interaction  problem.  By 
treating  the  oncoming  wakes  as  time-dependent  gusts,  these  models  attempt  to  model  how  the 
oncoming  wakes  periodically  change  the  incidence  of  an  isolated  airfoil.  None  of  these  models 
accounts  for  the  physics  of  the  chopped  wake  segments  within  the  blade  passage.  Meyer  [1958] 
and  Lefcort  [1965]  attempted  to  model  this  behavior.  They  modeled  the  chopped  wake  segments 
as  jets  sucking  fluid  from  the  suction  surface  of  a  compressor  blade  and  impinging  fluid  on  the 


13 


pressure  surface.  They  then  calculated  the  unsteady  static  pressure  on  an  isolated  thin  airfoil 
in  much  the  same  manner  as  Sears  [1941]  except  the  modeled  wake  segments  gave  large  local 
variations  in  static  pressure.  The  contribution  of  Lefcort  [1965]  allowed  for  a  more  realistic  wake 
width  in  the  model. 

So  far,  we  have  only  discussed  the  use  of  unsteady  thin  airfoil  theory  to  model  the  unsteady 
wake  interaction  problem.  Kemp  and  Sears  [1953]  and  more  recently  Lienhart  [1974]  have  used 
unsteady  thin  airfoil  theory  to  model  the  unsteady  potential  flow  interaction  as  well.  Kemp 
and  Sears  [1953]  attempted  to  solve  the  unsteady  interaction  problem  of  a  row  of  stator  blades 
followed  by  a  row  of  rotor  blades.  To  avoid  a  very  complex  mathematical  problem,  they  made 
several  assumptions.  The  unsteady  lift  of  a  single  stator  depended  only  on  the  relative  motion 
between  that  blade  and  the  steady  potential  flow  fields  of  the  rotor  blades.  The  unsteady  lift  of 
a  single  rotor  blade  depended  on  two  things.  First,  the  unsteady  lift  depended  on  the  relative 
motion  between  that  blade  and  the  steady  potential  flow  fields  of  the  stator  blades.  Also,  the 
unsteady  lift  on  each  rotor  blade  was  computed  using  the  upstream  stator  blade  wakes  with  a 
theory  similar  to  that  of  Sears  [1941]. 

Kemp  and  Sears  [1953]  developed  the  first  singularity  technique  for  use  with  multiple  blades. 
However,  when  calculating  the  unsteady  lift  of  a  particular  blade,  they  neglected  the  unsteady 
effects  from  all  other  blades.  As  the  solidity  of  a  blade  row  increases,  the  neighboring  blades  have 
a  stronger  influence  on  the  unsteady  lift  of  the  reference  blade.  We  can  make  a  slightly  better 
assumption  of  the  unsteady  lift  by  using  the  method  of  Weinig  [1935].  He  computed  the  values 
of  steady  lift  in  cascades  of  various  solidities  and  stagger  angles  which  allows  us  to  change  the 
value  of  (a  quasi-steady  correction).  This  method  does  break  down  for  high  stagger  angles, 
however. 

Some  researchers  have  developed  unsteady  cascade  analyses  using  unsteady  thin  airfoil  the¬ 
ory.  Whitehead  [  1960]  modeled  all  of  the  blades  with  a  finite  number  of  vorticies  (usually  eight) 
and  all  of  the  trailing  vortices  with  a  continuous  vortex  distribution.  Fleeter  [1973]  updated  this 
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method  to  include  the  effects  of  compressibility.  Whitehead  [1960]  also  considered  the  vibrational 
response  of  each  blade  to  the  unsteady  lift  and  Smith  [1971]  used  this  with  some  compressibility 
effects  to  determine  the  resulting  acoustic  field.  For  unsteady  cascade  analyses,  we  need  the 
inter-blade  frequency, 


The  strength  of  the  bound  and  trailing  vorticity  from  all  blades  must  be  summed  and  this 
sum  must  satisfy  the  unsteady  blade  boundary  condition.  In  the  development  of  his  analysis, 
Whitehead  [i960]  uses  a  summation  method  that  excludes  solutions  when  rn  equals  an  integer 
multiple  of  2tt.  Therefore,  his  analysis  does  not  solve  cases  where  the  number  of  cycles  in  the 
distortion  (the  inflow  harmonic)  is  a  multiple  of  the  number  of  blades.  After  some  complex 
integration,  Whitehead  [1960]  computes  an  unsteady  lift  coefficient  which  we  may  use  to  replace 
the  Sears  Function  in  the  unsteady  lift  equation  of  Naumann  and  Yeh  [1973]. 

Henderson  [1972]  also  developed  an  unsteady  cascade  analysis.  He  modeled  the  reference 
blade  with  a  continuous  vortex  distribution  and  the  other  blades  with  a  single  vortex  located 
at  the  quarter  chord.  Like  Whitehead  [1960],  Henderson  [1972]  also  modeled  the  trailing  vor¬ 
tices  with  continuous  vortex  distributions.  For  turbomachinery  problems,  we  should  note  that 
modeling  the  trailing  vorticity  along  the  extended  chord  line  can  become  a  poor  assumption.  Un¬ 
fortunately,  Henderson  [1972]  could  not  determine  an  analytical  solution  for  one  of  the  integrals 
in  his  analysis  and  therefore,  he  could  not  obtain  a  close-form  solution. 

Lienhart  [1974]  extended  the  potential  interaction  work  of  Kemp  and  Sears  [1953]  to  include 
some  vortex  representation  of  all  the  blades.  Thus,  when  calculating  the  unsteady  lift  on  a  blade, 
he  included  the  unsteady  effects  due  to  all  the  other  blades.  However,  he  did  not  include  unsteady 
wake  interaction.  Of  these  cascade  analyses,  only  Henderson  [1972]  includes  the  effects  of  camber 
or  steady  angle  of  attack. 

For  us  to  apply  the  preceding  singularity  techniques  to  the  turbomachinery  problem,  we  must 
calculate  the  unsteady  lift  at  several  blade  sections  in  two-dimensional,  independent  strips.  In  a 
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later  section,  we  shall  discuss  the  proper  technique  for  summing  unsteady  forces  and  moments 
in  a  turbomachine  using  two-dimensional,  unsteady  strip  theory. 

Three-dimensional  singularity  techniques,  although  mathematically  more  complex,  have  also 
been  developed.  Brown  [1964]  developed  an  unsteady  lifting  line  technique.  As  with  lifting  line 
methods  for  steady  flow,  a  single  vortex  models  each  blade  section  and  all  of  these  vorticies 
form  a  continuous  radial  vortex  distribution.  Again,  the  strength  of  this  distribution  changes  in 
time  with  the  varying  inflow  resulting  in  shed  vorticity  downstream  of  the  blade.  Brown  [1964] 
uses  corrections  to  account  for  the  bound  vortex  on  neighboring  blades  and  both  the  radial 
and  streamwise  trailing  vorticity.  For  high  reduced  frequencies  where  the  wavelength  of  the 
disturbance  is  much  smaller  than  the  chord,  we  might  suspect  that  modeling  each  blade  section 
with  a  single  vortex  will  not  give  as  enough  resolution  to  properly  calculate  the  unsteady  lift.  In 
order  to  somewhat  alleviate  this  problem,  we  could  model  each  blade  section  with  several  vorticies 
and  thus  model  each  three-dimensional  blade  with  several  lifting  lines.  We  call  this  technique  an 
unsteady  lifting  surface  technique.  Tsakonas,  Jacobs,  and  Ali  [1973],  Kerwin  and  Lee  [1978], and 
Tsakonas,  Jacobs,  and  Liao  [1983]  have  developed  such  unsteady  lifting  surface  techniques.  The 
latter  two  techniques  consider  moderately  loaded  blades  where  the  discrete  lifting  lines  lie  on 
the  camber  surface.  Kerwin  and  Lee  [1978]  also  include  sources  in  their  singularity  technique  to 
account  for  blade  thickness.  The  shed  trailing  vorticity  travels  on  a  helicoidal  path  downstream 
and  the  various  unsteady  lifting  surface  techniques  differ  on  how  to  model  this  complex,  three¬ 


dimensional  phenomenon. 
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Unsteady  Computational  Fluid  Dynamics 

Unsteady  singularity  techniques  have  several  drawbacks.  First,  we  cannot  adequately  model 
many  blades  of  interest  in  turbomachinery.  Singularities  cannot  give  us  good  surface  definition 
(unless  someone  develops  a  good  unsteady,  three-dimensional  panel  method)  and  most  methods 
limit  the  camber,  thickness,  and  steady  angle  of  attack.  We  also  have  difficulties  modeling  the 
shed  trailing  vorticity,  the  chopped  wake  segment  effects,  and  some  three-dimensional  effects. 
Finally,  since  the  singularities  can  only  induce  potential  flow,  we  cannot  extend  these  methods  to 
viscous  flow.  It  seems  that  much  of  our  future  unsteady  flow  computations  will  lie  in  numerically 
modeling  the  fluid  flow  instead  of  modeling  the  blades. 

Several  unsteady  numerical  methods  have  recently  been  developed.  After  much  initial  work 
by  several  researchers,  Hodson  [1985]  developed  a  numerical  method  to  compute  the  wake¬ 
generated  unsteady  flow  around  a  rotating  blade  row.  He  modeled  the  two-dimensional,  com¬ 
pressible,  Euler  equations  and  solved  them  using  an  explicit  finite  volume  approach.  Figure  6 
shows  some  typical  results  for  a  turbine  rotor  row.  These  results  show  much  of  the  anticipated 
physics  of  wake  interactions.  The  unsteady  static-pressure  distribution  may  be  numerically  in¬ 
tegrated  to  find  the  amplitude,  phase  angle,  and  location  of  the  unsteady  lift.  Koya  and  Kotake 
[1985]  developed  a  similar  method.  They  do  everything  the  same  as  Hodson  [1985]  except  that 
they  apply  the  method  to  the  three-dimensional,  compressible,  Euler  equations.  They  have  also 
attempted  to  apply  the  method  to  successive  blade  rows  to  obtain  the  unsteady  potential  flow 
effects.  However,  since  the  method  is  inviscid,  they  cannot  compute  the  wakes  from  the  upstream 
blade  row  that  interact  with  the  downstream  blade  row. 

We  still  have  many  problems  that  hinder  our  use  of  unste;.  \y  computational  fluid  dynamics  in 
computing  unsteady  forces  and  moments  in  turbomachines.  When  applying  Lulu  solvers  to  the 
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steady  cascade  problem,  we  use  spatial  periodicity  in  our  boundary  conditions.  In  the  unsteady 
cascade  problem,  we  can  either  use  time-dependent  periodicity  conditions  or  we  must  include 
the  number  of  blade  passages  over  which  the  upstream  flow  pattern  repeats.  These  alternatives 
are  either  difficult  to  apply  or  very  time  consuming.  Also,  we  must  generate  a  suitable  grid. 
For  the  multiple  blade  row  problem,  the  grids  for  each  blade  row  must  coincide  at  each  time 
step.  Additionally,  we  might  have  a  simple  bookkeeping  and  computer  storage  problem  since  the 
computations  must  include  enough  time  steps  to  include  the  transient  period  and  at  least  one 
cycle  of  the  periodic  unsteady  flow.  Finally,  no  adequate  computational  methods  exist  yet  for 
the  incompressible,  unsteady  flow  problem. 

Several  methods  exist  for  computing  unsteady  boundary  layer  flows.  Most  of  these  methods 
simply  use  the  unsteady,  two-dimensional,  boundary  layer  equations.  A  major  problem  occurs 
when  attempting  to  use  these  methods  in  turbomachinery  problems  since  these  boundary  layers 
have  regions  of  laminar,  transitional,  turbulent,  and  separated  flow.  We  find  the  biggest  problem 
in  our  selection  of  a  turbulence  model.  Modeling  Reynolds  stresses  is  difficult  enough  in  steady, 
viscous  flow.  For  unsteady  flow,  we  can  obtain  Reynolds  stresses  due  to  periodic  unsteadiness 
as  well  as  turbulent  unsteadiness.  Furthermore,  our  selection  of  empirical  constants  for  these 
turbulence  models  becomes  more  difficult  because  of  the  interaction  between  the  periodic  and 
turbulent  unsteadiness.  Besides,  few  if  any  good  periodic  unsteady  experiments  exist  for  the 


selection  of  these  constants. 
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Nonuniform  Inflow 

Before  we  compute  the  unsteady  forces  and  moments  within  a  turbomachine,  we  need  to 
establish  the  nonuniform  flow  upstream  of  the  blade  row  of  interest.  We  can  compute  the 
nonuniform  inflow  by  numerically  solving  the  viscous  flow  around  the  blades  creating  the  wakes, 
or  we  can  use  some  numerical  method  in  conjunction  with  some  empirical  wake  models.  However, 
in  order  to  satisfactorily  obtain  the  entire  harmonic  content  of  the  nonuniform  inflow,  we  may 
need  to  measure  the  wakes  experimentally.  In  any  case,  we  must  be  aware  that  the  flow  field 
experienced  by  a  rotating  blade  row  will  be  different  than  the  flow  field  that  exists  with  the 
rotating  blade  row  removed.  Since  the  rotor  blades  accelerate  the  flow,  the  rotor  blades  will  shift 
the  streamlines  toward  the  roots  of  the  rotor  blades  'or  to  satisfy  continuity. 

If  we  intend  to  compute  the  uni  «ady  forces  and  moments  with  a  singularity  technique,  we 
need  to  simulate  the  nonunuo  m  inflow  as  a  convected  gust.  We  can  consider  any  nonuniform 
inflow  as  a  harmonic  function  in  space.  First,  we  decompose  the  nonuniform  inflow  into  a 
circumferential  average  velocity  and  a  circumferentially  varying  velocity, 

t'(z,r,0)  =  V(z,r)  +  V(z,r,0)  . 

We  may  describe  this  velocity  by  a  Fourier  series  as 

OO 

V(z,r,  0)  =  ^  ^  [an  cos(n^)  +  bn  sin(n0)] 

n  =  l 

where 

an(z,r)  =  —  f  V  cos (n9)d9 
b„(z,r)  =  —  f  V  sin(n0)d0  . 

*  J-w 


and 
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The  term  represents  the  circumferential  average  velocity,  V .  We  can  now  rewrite  the  distortion 
velocity  for  the  nth  harmonic  as 

Vn(z,r,0)  =  an  cos (nO)  +  bn  sin(n0) 

=  \Vn\{a„/\Vn\)  cos(n9)  +  \Vn\ (bn/\Vn\)  sin(n0) 

=  \Vn\  cos{<t>n)  cos(n0)  +  \Vn\ sin(<£„)  sin(nfl) 

=  |K, |  cos(n0  -  4>n) 

where  |  V„  |  represents  the  amplitude  of  the  distortion  velocity  and  <f>n  represents  the  phase  angle. 
The  phase  angle  represents  the  combination  of  cosine  and  sine  waves  that  make  up  the  nth 
harmonic.  We  can  define  the  amplitude  and  phase  angle  as 

|Vnl(*,r)  =  sjal  +  bi 


and 

<f>n(z,r)  =  tan_1(6n/a„)  . 

We  define  the  phase  angle  as  positive  when  measured  in  the  direction  opposite  to  blade  rotation. 
It  is  also  important  to  know  the  6  position  where  the  measured  phase  angle  was  referenced.  A 
typical  reference  position  might  be  vertically  up.  We  must  also  remember  that  all  amplitudes 
and  phase  angles  can  vary  radially  and  these  variations  can  lead  to  changes  in  unsteady  forces. 
For  simplicity,  we  will  rewrite  the  distortion  as  an  exponential, 

V(z,r,9)  =  |Kn|ei(n*  "  *»>  . 


In  writing  the  distortion  in  this  manner,  we  are  only  concerned  with  the  real  part  of  the  expression. 
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Unsteady  Forces  and  Moments  on  a  Rotating 
Shaft 


We  will  now  focus  on  the  calculation  of  unsteady  forces  and  moments  on  a  rotating  shaft 
caused  by  rotor  blades  interacting  with  a  nonuniform  inflow.  In  order  to  accomplish  this  cal¬ 
culation,  we  will  need  to  transfer  the  nonuniform  inflow  into  a  rotating  coordinate  system.  For 
the  interaction  of  a  stator  blade  row  with  upstream  rotor  wakes,  the  rotating  rotor  wakes  can  be 
transferred  into  a  stationary  coordinate  system.  In  either  case,  we  need  to  adjust  the  distortion 
velocity  in  order  to  apply  a  two-dimensional,  unsteady  thin  airfoil  theory  to  a  blade. 


The  distortion  velocity  will  depend  on  the  meridional  direction1,  m.  The  fluid  flows  merid¬ 
ionally  a  distance  m  in  time  In  this  time,  the  fluid  also  moves  circumferentially  a  distance 

'  m 

rAO.  Moving  rAO  takes  time  5^.  Since  these  two  times  must  be  equal,  we  obtain 


A  e 


mV  e 

'rVZ 


Thus,  if  swirl  exists  upstream  of  the  rotating  blades,  the  phase  angle  will  vary  meridionally: 

_  »  fn  (t  -  =£*■)  -  *„] 

Vn(z,r,6)  =  \Vn\e  l  V  J  . 


We  can  estimate  the  m  dependence  of  |K„|  with  a  factor  e~6nmlVm.  6n  is  a  factor,  in  units  of 
l/sec,  that  shows  how  the  distortion  amplitude  changes  as  m  increases.  One  example  of  how 
|l/„|  can  change  meridionally  is  the  diffusion  of  a  wake  due  to  viscosity.  Thus 


Vn(z,r,0)  = 


IK,  | 


0-6nm/V 


v  <(n9' 


/ 


1  The  meridional  direction  is  defined  by  the  direction  of  the  vector  sum  of  the  axial  and  radial 
mean  velocities. 
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Since  most  analytical  techniques  compute  the  unsteady  lift  on  an  isolated  airfoil,  or  at  least 
reference  the  computation  to  a  particular  blade,  we  shall  include  a  shift  in  the  0  position  due  to 
the  rotor  blade  number,  B, 


Vn =  \Vn\e~6'mlv ~  e 


y  *  n(9  +  2Si) 

"e  1 


n  mV g 

rV  m 


In  order  to  calculate  the  unsteady  lift  on  a  rotor  blade,  we  must  shift  from  the  stationary 
coordinate  system  that  we  have  been  using  to  a  coordinate  system  rotating  with  the  blades. 
Using  Figure  7,  we  see  that  we  must  transform  from  ( z ,  r,  6)  to  (z',  r',  6')  by 

e  =  0'  +  tl>mc  -  n t, 

z  =  z>  +  2rmc  , 

and 

r  =  r  . 

ip  me  is  the  skew  angle  of  the  blade  midchord,  Q  is  the  angular  velocity  of  the  blades,  and  zTmc  is 
the  rake  of  the  blade  midchord.  The  reference  point  of  Figure  7  could  represent  a  radial  stacking 
line  or  a  radial  line  defining  the  midchord  rake  and  skew  to  be  zero  at  the  hub.  We  define  the 
skew  angle  as  positive  in  the  direction  opposite  to  blade  rotation  and  rake  as  positive  in  the  same 
direction  as  axial  flow.  Assuming  that  zr„c  is  in  the  meridional  direction  instead  of  the  axial 
direction  (since  zrme  is  usually  small),  we  can  now  write  the  distortion  velocity  in  terms  of  the 
rotating  coordinate  system 

Vn(z',r',e',i,£)  =  |Vn|e-Mm' +  ^mc)/Fm 

cf(-"n‘  +•&-*.  +  -  70) 

In  order  for  us  to  view  the  problem  as  time  increases,  we  can  simply  change  the  sign  of  all 
the  terms  in  the  last  exponential  function.  This  will  not  change  our  results  since  we  are  only 
concerned  with  the  real  part  of  the  expression.  This  yields 
Vn(z'y,0',t,t)  =  \Vn\e-s^m'  + 

-  ni&mc  +  "  -  »*'  +  v^) 


We  note  here  that  the  measured  phase  angle,  4>n,  causes  a  time  lead  and  the  skew  angle  causes 


a  time  lag  in  the  distortion  relative  to  the  case  with  no  measured  phase  and  no  skew.  The  rake, 
however,  causes  a  time  lead  if  V»  is  positive  and  a  time  lag  if  Vg  is  negative. 

Many  analytical  methods  for  computing  the  unsteady  lift  use  a  coordinate  system  (x'b,  y'b,  r') 
attached  to  the  rotor  blade  as  seen  in  Figure  7.  We  can  rotate  our  current  coordinate  system  to 
this  rotating  system  using  the  following  relations: 

m!  =  x[  cos(7  +  o)  -  t/b  sin(7  +  a) 


and 

r'0'  =  x'b  sin(7  +  a)  +  yb  cos(7  4-  a)  . 


From  the  velocity  diagram  of  Figure  7,  we  obtain 


sin(7  +  a) 


Vb  +  r'Q 

W 


and 


and  these  equations  yield 


cos(7  +  a) 


m 


x[V_m  _  i/t  (VeJ  r'Q) 
W  W 


and 


=  (V*  +_  r'n)  yjVm 

r'W  r'W 

From  the  velocity  diagram  of  Figure  7,  we  get 


Ve 

V _ 


tan(/?)  . 
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We  can  now  rewrite  the  last  two  terms  in  the  last  exponential  of  the  expression  for  the  distortion 
velocity, 


_n0'  +  nm'v<>  _  nxb(V»  +  r'ty  _  ny>bvm  nx'bVe 

"  r'Vm  r'W  r'W  r'W 

__  ni/bVe(ye+'  r'ft) 
r'VmW 


ndx'b 

W 


nClxb 

W 


W 

W 


Ve(V0  +  r'Q)  +  vmv„ 
r'fi 

VmVm  +  Ve(ye  +  r'fl) 
-VmVe  +  r'QVm  +  VmVe 


nfiarj  nfij/j 

'l  +  I 

(fc)l 

'VtJz  r'n\ 

i  V  m  / 

W  W 

-1 

(Zi.) 

\vm) 

1  + 

(2^0) 

nflxj  nflj/j  I"  1  —  tan  (/?)  tan  (7  + 
~  W  W  L  tan  (/ 3 )  +  tan  (7  + 

nflxl  nO  j4  ...  . 

= - =2-  —  — *  cot  (/?  -f  7  -f  a) 

W  W  v^  7  ' 

/u\  nfVb 

w  v*v  w 


•<»)' 
«)  . 


We  now  obtain 

?«(*»,  1/1, r\*,*)  =  \Vn\e-s^m'  + 


;["««  -***  +  ♦.-  +  (f)^] 


The  calculations  of  the  unsteady  lift  in  the  (x'b,  \/b,  r')  coordinate  system  requires  us  to  resolve 
the  distortion  velocity  into  components  in  the  x'b  direction,  u,  and  the  yj,  direction,  v.  We  can 
compute  these  components  as  follows 

un(x'b,yb,r',t,e)  =  V„  cos(/3  +  7  +  a) 

=  Vn  cos(/3)  cos(7  +  a)  -  V„  sin (0)  sin(7  +  a) 

~  ccs(7  -f  a)  I  sir. (7  +  a) 

and  _ 

Vn(x'h,y'b,r' ,t,£)  =  -V„  sin(/3  +  7  +  o) 

=  -Vn  sin(/?)  cos (7  +  a)  -  lr„  cos(/7)  sin(7  +  a) 

=  -Vmn  sin(7  +  a)  +  cos(7  +  o)  . 
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In  this  section,  we  will  use  the  equation  for  the  unsteady  lift  per  unit  span  on  an  isolated 
airfoil  derived  by  Naumann  and  Yeh  [1973], 

Ln  =  )~pcw\vnS(un)  +  aunFa(un)  +  F}(wA  . 

i  act  c  J 

The  expressions  for  u  and  v  were  used  in  the  derivation  of  this  equation.  The  term  involving  ijb 
dropped  out  since  we  are  only  interested  in  the  effects  on  the  blade  surface  where  x/b  =  0  in  the 
thin  airfoil  approximations.  The  term  involving  x'b  became  part  of  three  larger  expressions  in 
deriving  the  functions  S,  Fa,  and  Fj .  These  expressions  were  integrated  from  x'b  =  —  |  to  x'b  =  |. 
For  the  chord,  we  should  estimate  the  projection  of  the  chord  in  the  meridional  direction.  This 
allows  us  to  use  a  two-dimensional  method  along  a  circumferential-mean  streamsurface  where 
three-dimensional  flow  effects  should  be  minimised.  The  values  of  c,  W,  a,  and  ym  are  all 
a  function  of  the  radius. 

Before  computing  the  total  unsteady  force  components  due  to  all  blades,  all  radial  sections, 
and  all  inflow  harmonics,  we  must  decide  where  the  unsteady  lift  acts  on  each  blade  section.  In  his 
analysis,  Sears  [1941]  calculated  that  the  unsteady  lift  would  act  at  x'b  —  —  |.  We  must  calculate 
the  unsteady  lift  at  the  location  in  which  it  acts  in  order  to  avoid  calculating  the  corresponding 
unsteady  moments  along  the  chord.  For  small  values  of  a,  the  relationship  between  the  midchord 
skew  and  rake  and  the  skew  and  rake  at  x'bo  is 

x'b 

Vw  «  ^  ~  —T  sin(7) 
r 

and 

Zmc  «  ~  *4.  C0S(7)  . 
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Using  Figure  7  and  the  law  of  sines  gives 

W  sin(f  -/?)  _  cos(/g) 

r'ft  sin(/?  +  7  +  a)  sin(/?  +  7  +  a) 

The  unsteady  thin  airfoil  theory  that  we  use  to  compute  the  unsteady  lift  assumes  a  small  value 
of  a.  Using  a  small  angle  approximation  gives 

sin(/?  +  7  +  a)  fts  sin(/?  +  7)  =  sin(/?)  cos(7)  +  cos(/?)  sin(7)  . 

FVom  the  velocity  diagram  of  Figure  7,  we  get 

~w~  ~  tan^ 

'  m 


and  finally 


+  W 


,  ,  nzrVg 
-nip  +  — = — 
r'V  m 


+ 


f 

W  \ 


cos(P)  sin(7) 


W  1  sin(/?)  cos(7)  +  cos(/?)  sin(7) 


—  sin(/?) 

cos(/?)  cos(7)  ( 

cos(/3) 

sin(/?)  cos(7)  +  cos(/?)  sin(7)  J 

.  nzrV  g  nSlx'b 
«  -ni/>  +  _  +  . 

r'Vm  W 

This  equation  yields  two  methods  to  compute  the  unsteady  lift  at  the  quarter  chord.  We  can  move 
the  coordinate  system  from  (x'b,  t/j,  r')  to  (x'b  —  |,j /b,r')  and  simply  use  the  values  of  skew  and 
rake  at  this  new  location.  Or  we  can  compute  the  unsteady  lift  at  midchord  and  then  move  the 
location  of  the  unsteady  lift  to  the  quarter  chord  by  multiplying  by  e-,rtnc/4W .  The  dependence 
of  the  location  of  the  coordinate  system  in  computing  unsteady  lift  can  be  demonstrated  by  the 
lift  function  developed  by  Brown  (1964]  as  seen  in  Figure  8.  For  a  particular  reduced  frequency, 
the  amplitude  of  unsteady  lift  does  not  change  with  coordinate  system  location.  Only  the  phase 
angle  of  the  unsteady  lift  is  affected.  We  can  compute  the  change  in  unsteady  lift  by 


=  Ln 


infix  i/W 
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Since  we  must  eventually  sum  the  unsteady  forces  radially,  we  should  compute  the  unsteady 
lift  at  the  location  where  the  unsteady  lift  acts.  This  new  rotating  coordinate  system  is  shown 
in  Figure  9  and  the  transformation  from  the  stationary  coordinate  system  to  this  new  rotating 
coordinate  system  would  have  used  the  skew  and  the  rake  at  quarter  chord  instead  of  at  midchord. 
Therefore,  to  obtain  the  unsteady  lift  acting  at  quarter  chord,  we  use  the  following  expression 
for  the  distortion  velocity  prior  to  using  the  equation  for  the  unsteady  lift: 

Vn(x'b,y'b,r',t,e)  =  + 

e‘h  #  *  (o^j 

After  adjusting  the  measured  inflow  phase  angle  for  skew  and  rake  at  quarter  chord,  we 
compute  an  unsteady  lift  for  each  blade, 

Zn(r\f,*)  =  |Z„|e*(nn‘ 

This  unsteady  lift  acts  in  the  positive  y*b  direction  (see  Figure  9).  Now  we  can  determine  the 
components  of  this  unsteady  lift  in  the  z' ,  r',  and  6'  directions  as 

Fz'n(r' ,t,£)  =  -L„  sin(7  +  a)  cos(/i)  , 

Fr‘n(r',t,£)  =  -Z„  sin(y  +  a)  sin(/i)  , 

and 

Frn{r',t,Z)  =  Z„cos(7  +  a) 

where 

“ =  un"  (£)  - 

fi  represents  the  inclination  of  the  streamlines  to  the  axial  direction.  To  calculate  the  total 
unsteady  force  acting  on  the  shaft  in  the  z'  direction,  we  must  move  all  of  the  Fz<n  vectors  to  a 
point  on  the  shaft  centerline.  This  involves  summing  the  Ft vectors  over  all  inflow  harmonics, 
all  blades,  and  all  radial  sections, 

mo  =  - 

r1  =  r h  n  =  l 
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FYom  the  summation  relationships  of  Gradshteyn  and  Ryzhik  [1980],  we  can  easily  derive 


( B  if  n  =  kB 
1 0  if  n  ^  kB 


Using  this  relationship  allows  us  to  obtain 


r-p  oc 

Fz'{t)  =  -B  ^2  l^tsf  sin(7  +  a)  cos(fx)  e,(kBnt  +  **s)  . 

r'=ir h  k  =  \ 

This  force  is  the  negative  of  the  unsteady  thrust  on  the  shaft.  The  unsteady  thrust  occurs  at 
frequencies  that  are  integer  multiples  of  the  blade  passing  frequency  and  can  be  computed  from 
the  inflow  harmonics  corresponding  to  these  frequencies.  We  denote  the  multiples  of  blade  rate 
harmonics  with  the  subscript  k. 


We  can  also  calculate  the  unsteady  side  forces  on  the  shaft.  The  components  of  unsteady 
force  in  the  r'  and  O'  directions  must  be  resolved  into  components  of  unsteady  force  in  the  x' 
and  y1  directions.  Using  Figure  10,  we  calculate  the  total  unsteady  force  in  the  x'  direction, 

mo  =  E  EE  |Z„|sin(7  +  <*)  sin(/*)  cos  (if)  -  ~ 

r'=r«i=lnrl  ' 

e.(«m  -  +  *„) 

tt  B  oo 

-  E  EE  \Ln\ cos(7  +  o)  cos 

r'=rH  Isl  n=l 
ei(ntU  -  +  *.) 


We  can  write  the  last  cosine  factor  in  each  term  as  a  sum  of  exponentials  using  Euler’s  Formula, 


^  53  ££Msin(7  +  «)sin(,) 

r> sr«  /=1  n=l 

e<[r»n<  -  (n+l)^  +  4-„  +  -  f] 

.  tt  B  oo 

„|sin(7  +  a)  sin(//) 

r'=rn  t=  1  f»=l 

e<[nn«  -  (n-l)2^  +  4>„  -  ^  +  f] 

1  TT  B  00 

-  2  53  5353  l^n|cos(7  +  a) 

r'—rj-i  /=rl  rt= 1 

e«[nn<  -  (n+l)2^  +  +  <£] 

.  rr  B  oo 

-  2  53  5353iz"icos(^+a) 

r'=rH  /=1  r>=ol 

e.(nn<  -  +  *„  -  0]  _ 

The  summation  relation of  Gradshteyn  and  Ryihik  [1980]  yield 

V''  -»(n±i)2*7/B  _  f  &  **  (nil)  =  kB 
jr[  \0  if  (nil)  £  kB 

We  can  now  write 

R  rr  oo 

=  7  5  53l^iB-1l  [«“(7  +  «)  sin(/i) e-,>/2  -  cos(7  +  a)] 

r'  zzr h  k  —  l 

ei[(*B-l)nt  +  QkB-i  +  V1] 

B  rr  °°  ~  [  1 

+  y  53  53|Z*a+i|[sin(7  +  «)  s^)«‘r/2  ~  cos(7  +  a)J 

r'—rn  k  =  l 

eH(kB+l)(U  +  4i.b+i  -  V>] 

We  can  compute  the  total  unsteady  force  in  the  t/  direction  in  the  same  manner: 

Fy'(t)  =  Y  53  53  1  -  sin(7  +  O')  Sin(/i)  -  cos(7  +  a)  e-'^2] 

r'=rH  k=l 

ei[(kB-l)nt  +  *tB-i  +  tf] 

+  j  53  53}£*3+i!  [-sin(-y  +  «>  sin(^)  -  cos(7  +  a)  e,,r/2] 
r1  =zth  k  =  1 

c«[(fcB+i)n<  +  4>fcfl+i  - 
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In  order  to  better  compare  Fy<  with  Fxi,  we  can  also  write  Fy>  as 

R  fT  00 

Fy'(t)  =  j  [sin(7  +  a)  sin(p)  e_,,r/2  -  cos(7  +  a)] 

r‘—rn  k  =  l 

gi[(fc£  —  l)Qt  +  +  ^  —  -j] 

B  rT  00 

+  J  £lWl[sin(7  +  «)  sm(/i)eW2  -  cos(7  +  a)] 

t‘=th  i  =  l 

ei[(itB+l)nt  +  *ke+i  -  <P  +  f ] 

In  the  rotating  coordinate  system,  the  unsteady  side  forces  occur  at  frequencies  adjacent  to  the 
integer  multiples  of  blade  passing  frequency  and  are  computed  from  inflow  harmonics  correspond¬ 
ing  to  these  frequencies.  A  90  degree  phase  shift  exists  between  the  portions  of  unsteady  side 
force  in  the  x'  and  r/  directions  that  occur  at  each  of  the  adjacent  frequencies. 

The  unsteady  forces  on  the  blades  also  create  unsteady  moments  at  the  point  on  the  shaft 
centerline  where  we  computed  the  unsteady  forces.  Vectorially,  we  can  compute  the  unsteady 
moment  at  this  point  by 

M(r',t,e )  =  x  Ln(r',t,£)  . 

The  components  of  this  unsteady  moment  in  the  z',  r',  and  O'  directions  are 

=  r'FK 

=  r'Ln  cos(7  4-  a)  , 

=  -zrFe'n 


=  r' Ln  sin(7  +  a)  cos (fi)  —  zrLn  sin(7  +  a)  sin(/i)  . 

We  can  compute  the  components  of  the  total  unsteady  moment  in  the  z',  x',  and  j (  direction  in 

the  same  manner  as  the  unsteady  force  components: 


AMO  =  B  £  jr\ZkB\r'  cos(7  +  a) +  , 

t'=th  *  =  1 
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_  R  T  ~  r 

MX‘(t)  =  -  ^  ^|Z*B-il  Ur  cos(7 +  a)e_,,r/2  -  r'  sin(7  +  a)  cos(^) 

r'=Tn  (fc=l 

-f  zr  sin(7  +  a)  sin(/x)j  e'[(*B-i)m  +  1  +  ^1 

B  rT  00  ~  r 

+  Y  cos(7  +  a)e,,r/2  -  r'  sin(7  +  a)  cos(fi) 

r'=rn  k'=  1 

+  rr  sin(7  +  a)  sin(/i)]  ^K^+Uni  +  ***+>  ~  *]  , 


* _  B  ~  r 

My/(<)  =  —  Y2  cos(7  +  a)e-/2  _  r' sin(7  +  a)  cos(/i) 

r'  =rH  k  =  l 

+  zr  sin(7  +  a)  sin(p)]  +  *‘B->  +  *  ~  f] 

B  rT  00  ~  r 

+  T  H  X^*B+1lrr  cos(7  +  <*)e,,r/2  -  r' sin(7  +  a)  cos(/i) 
r'=rH  *  =  1 

+  rr  sin(7  +  a)  sin(//)j  e,'[(lB+1)n‘  +  *‘B+>  "  *  +  f ]  . 

The  unsteady  torque,  — Af,<,  occurs  at  integer  multiples  of  the  blade  passing  frequency  while 
the  unsteady  bending  moments  occur  at  the  frequencies  adjacent  to  the  integer  multiples  of  the 
blade  passing  frequency. 

We  have  now  computed  the  various  components  of  unsteady  force  and  unsteady  moment  on 
the  rotor  shaft  for  a  coordinate  system  rotating  with  the  shaft.  Next,  we  need  to  transform  these 
forces  and  moments  to  a  stationary  coordinate  system  as  shown  in  Figure  11.  The  z  components 
of  the  unsteady  force  and  moment  are  exactly  the  same  as  the  z'  components: 

Fz(t)  =  -B  f^f2\LkB\sm('(  +  a)cos(fi)e^kBnt  +  *kB'> 
r=rH  k  =  1 


Mz(t)  =  B  ^|Zfcflir  cos(7  +  a)e,(tsn<  +  *‘s)  . 


r= rH  k  =  1 


Using  Figure  12,  we  can  write  the  relation  between  the  two  coordinate  systems  as 


x  =  x'  cos(— Q/)  +  y'  sin(— fi<) 


y  =  —  x'  sin(— Qt)  -f  y'  cos (— • 
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The  transformations  for  unsteady  side  forces  are 

Fx(t)  =  Fx>  cos(—  fit)  4-  Fy>  cos 


=  Fxi  cos(fii)  4-  Fyi  cos 


Fy(t)  =  —  Fxi  cos  4-  cos(-fii) 

=  -Fr/  cos  4-  4-  Fy<  cos(Q<)  . 


We  can  change  these  transformations  to  exponentials  using  Euler’s  Formula: 

Fx(t)  =  i[Fr-e*'n<  +  Fx‘  e~iClt  +  Fy,e’(nt+S)  4-  Fy  e-,(n‘  +  *)] 


#V(f)  =  i  [-F#.ei(n,+  *)  -  Fx-e-(n‘+  f)  +  Fy,e<n‘  4-  Fy- e— nt]  . 

We  can  now  compute  the  unsteady  horizontal  side  force  in  the  stationary  coordinate  system: 

rT  °o 

hw  =  T  £ 


r=rH  *  =  1 


->'*/ 2  __ 


^sin(7  +  a)  sin(/z)e-,,r/2  -  cos(7  4-  a)|  e,iBn‘ 

+  [sin(7  +  a)  sin(/i)e“”r/2  -  cos(7  4-  a)j  e,(iB~2)n< 

4-  [sin(7  +  a)  sin(/i)  e-”^2  —  cos(7  4- »)]  e,kBClt 
4-  J-sin(7  +  a)  sin(/i)e-,,r^2  4-  cos(7  +  a)j  e‘(tB-2)n<  | 


p«(*kB-l  +  VO 


+  T  £  SlZ^+i 


r=rH  k=  1 


,iir/2  _ 


<(fcB+2)at 


jsin(7  +  a)  sin(ji)  e”/2  -  cos(7  4-  a)j  e‘^B+2^n< 

4-  |sin(7  4- a)  sin^)6'^2  —  cos(7  4-  a)j  e'kBn> 

4-  J-sin(7  4- a)  sin(^)e,,r/2  4-  cos(7  4-  a)j  e,(tB+2)n< 
4-  sin(7  4- a)  sin(/i)e,ir^2  -  cos(7  4-  a)J  etfcBrit  \ 


p»(*kB  +  l  -  VO 
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This  expression  simplifies  to 

~  B  TT  00  ~ 

Fx(t)  =  —  £  £  l^ts-il  [sin(7  +  «)  sin(fi)  e~i,r/2  -  cos(7  +  <*)] 

r=rjj  k= 1 

ei(kBClt  +  4>*b_i  +  <//) 
g  rT  00 

+  2  E  £|£*b+iI  [sin(7+a)sin(/i)e^2  -  cos(7  +  a)] 

r  =  rH  jfc  =  i 

ei(kBSU  +  4>kB  +  i  - 

We  can  formulate  the  expressions  for  Fy,  Mz,  and  My  in  the  same  manner: 

B  TT  00  ~  r 

^v(0  =  y  £  £lW.il[sin(7  +  «)  sin(^)c-'/2  -  cos(7  +  a)j 

r=r//  fc  =  l 

ei(kB(it  +  tke-i  +  V>  -  f) 

5  rT  00  ~  r 

+  Y  £  £|£*fl+il  [sin(7  +  a)  sin(p)e‘^2  -  cos(7  +  a)l 

r=rH  k  =  l 

ei(kBiit  +  tjs+i  -  <li  +  -f) 

_  B  rT  ~  r 

Mx(t)  -=  y  £  £|ZjtB-i|  Ur  cos(7  +  a)  e_,,r/2  -  r  sinfr  +  a)  cos(^) 

r=rH  *  =  1 

+  zr  sin(7  +  a)  sin(/i)]  +  *«-*  + 

B  rT  00  ~  r 

+  y  £  £|ZjfcB+i|[^r  cos(7  +  a)e,,r/2  -  r  sin(7  +  a)  cos(/x) 

r=rH  k  =  l 

+  zT  sin(7  +  a)  sin(p)J  e'^kBnt  +  **B+>  _  ^  , 

and 

A fy(t)  =  —  £  £|Z*b_i|  zr  cos(7  +  a)e-,,r/2  -  r  sin(7  +  a)  cos(/i) 

r=rH  *  =  1 

+  zr  sin(7  +  a)  sin(/i)]  e<kBnt  +  *‘B-  +  *  ~  i) 

B  rT  00  r 

+  T  £  £|Z*B+ilk  cos(7  +  a)e"'2  -  r  sin(7  +  a)  cos(fi) 

t=th  k  =  l 

+  zr  sin(7  +  a)  sin(/i)J  e’(kBn‘  +  *iB+‘  -  ^  +  f )  . 

The  unsteady  side  forces  and  bending  moments  are  exactly  the  same  in  the  rotating  and  station¬ 
ary  coord'uate  systems  except  that  they  occur  at  frequencies  adjacent  to  the  integer  multiples  of 
blade  passing  frequency  in  the  rotating  coordinate  system  and  they  occur  at  the  integer  multiples 
of  blade  passing  frequency  in  the  stationary  coordinate  system. 
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Conclusions 


We  have  presented  a  method  for  computing  the  unsteady  forces  and  moments  on  the  shaft  of 
a  rotor  blade  row  operating  within  a  nonuniform  inflow.  The  method  uses  a  Fourier  analysis  and 
several  coordinate  transformations  to  develop  equations  for  all  six  components  of  unsteady  forces 
and  moments-in  both  a  stationary  coordinate  system  and  a  coordinate  system  rotating  with  the 
rotor  shaft.  This  strip  theory  does  incorporate  the  rake  and  skew  of  the  rotor  blades.  In  the 
stationary  coordinate  system,  all  of  the  components  of  the  unsteady  forces  and  moments  occur  at 
integer  multiples  of  the  blade  passing  frequency.  The  unsteady  thrust  and  torque  are  computed 
from  the  inflow  harmonics  corresponding  to  these  frequencies  while  the  unsteady  side  forces  and 
bending  moments  are  computed  from  the  inflow  harmonics  corresponding  to  frequencies  adjacent 
to  the  integer  multiples  of  blade  passing  frequency.  The  only  difference  in  these  results  when 
we  switch  to  a  coordinate  system  rotating  with  the  blades  is  that  the  unsteady  side  forces  and 
bending  moments  actually  occur  at  frequencies  adjacent  to  the  integer  multiples  of  blade  passing 
frequency. 

This  two-dimensional  strip  theory  can  be  used  in  conjunction  with  any  two-dimensional 
method  for  computing  the  unsteady  lift  on  a  single  rotor  blade  operating  within  a  nonuniform 
inflow.  If  we  use  an  unsteady  thin  airfoil  technique,  the  method  proves  to  be  very  fast  and 
adequate  for  predicting  trends.  A  computer  program  has  been  written  using  the  unsteady  thin 
airfoil  theory  of  Naumann  and  Yeh  [1973]  with  options  to  use  the  methods  of  Weinig  [1935]  and 
Whitehead  [I960].  This  Penn  State  ARL  code  uses  Fourier  analysed  wake  surveys  measured 
at  several  radii.  Because  of  the  finite  number  of  radial  and  circumferential  measurements  as 
well  as  the  many  approximations  involved  in  the  unsteady  thin  airfoil  technique,  we  should 
not  depend  on  these  calculations  to  obtain  absolute  unsteady  force  levels.  However,  comparing 
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calculations  using  different  blade  designs  gives  important  trends  of  how  the  different  designs 
affect  the  unsteady  forces  and  moments.  Hopefully,  engineers  will  incorporate  these  and  other 
unsteady  force  computations  within  the  turbomachinery  design  process. 
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Figure  I.  Periodic-average  wake  transport  and  interaction. 
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Figure  3.  Transverse  distortion  velocity  (Sears [ 1941 ]) . 
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Figure  6.  Wake-generated  unsteady  flow  around  a  rotating  turbine 
blade  row. 
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Figure  10.  System  used  to  resolve  the  unsteady  lift  force  into 
components  in  the  x '  and  y ’  directions  (flow  out  of 


Figure  11.  The  stationary  coordinate  system. 
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